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The concordance of standard big bang nucleosynthesis theory and the related observations of the light element 
isotopes (including some new higher *He abundances) will be reviewed. Implications of BBN on chemical evolution, 
dark matter and constraints on particle properties will be discussed. 
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1. Introduction 

The question of the concordance of standard 
big bang nucleosynthesis (BBN) rests on our abil- 
ity to compare the predictions of the standard 
model with the available observational data. At 
this time, as I will argue below, there is indeed 
concordance between the predictions of BBN and 
the observational determination of the primordial 
abundances of ^He, ''Li, and some measurements 
of D. Our lack of understanding of the chemical 
and stellar history of ^He prevents us from using 
this isotope as a critical check of BBN. Recent 
measurements of the '''He mass fraction which in- 
dicate a slightly higher abundance, diminish the 
role of ^He as a critical test as the higher abun- 
dances are very insensitive to the value of the 
baryon-to-photon ratio, t], the key undetermined 
parameter of BBN. In this case, we must wait 
until a definitively reliable abundance of D/H (if 
possible) is available from quasar absorption sys- 
tems to pin down the value of rj and definitively 
check the concordance of BBN. 

At this time, it is nevertheless still possible to 
draw some important inferences from BBN con- 
cerning the chemical evolution of the light ele- 
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ments, the amount of non-baryonic dark mat- 
ter, as well as constraints on particle properties. 
These issues will be touched upon briefly here. 
For a more complete discussion see ||^. 

BBN took place in the early universe when the 
temperature scale was < 1 McV. The key events 
leading to the synthesis of the light elements fol- 
lowed from the period when the weak interaction 
rates were in equilibrium, thus fixing the ratio 
of number densities of neutrons to protons. At 
T ^ 1 MeV, (n/p) ~ 1. As the temperature 
fell and approached the point where the weak 
interaction rates were no longer fast enough to 
maintain equilibrium, the neutron to proton ratio 
was given approximately by the Boltzmann fac- 
tor, (n/p) ~ e"'^™/"^, where Am is the neutron- 
proton mass difference. The final abundance of 
^He is very sensitive to the (n/p) ratio. 

The nucleosynthesis chain begins with the for- 
mation of deuterium through the process, p+n — > 
D +"f. However, because the large number of pho- 
tons relative to nucleons, i]^^ = n~^/nB ^ 10^°, 
deuterium production is delayed past the point 
where the temperature has fallen below the deu- 
terium binding energy, Eb — 2.2 MeV (the aver- 
age photon energy in a blackbody is Ej ~ 2. IT). 
When the quantity r]^^exp{—EB /T) ^ 1, the rate 
for deuterium destruction (D +7 —tp + n) fi- 
nally falls below the deuterium production rate 
and the nuclear chain begins at a temperature 
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T ~ O.lMeV. 

The dominant product of big bang nucleosyn- 
thesis is ^He, resulting in an abundance of close to 
25% by mass. In fact, it is of key importance that 
there are no observations of any any system for 
which the ^He mass fraction lies significantly be- 
low 25% thus indicating a primordial origin for 
this element. In the standard model, the ''He 
mass fraction depends primarily on the baryon to 
photon ratio, rj as it is this quantity which deter- 
mines the onset of nucleosynthesis via deuterium 
production. But because the (n/p) ratio is only 
weakly dependent on r], the ^He mass fraction is 
relatively flat as a function of rj. When we go 
beyond the standard model, the ''He abundance 
is very sensitive to changes in the expansion rate 
which can be related to the effective number of 
neutrino flavors as will be discussed below. Lesser 
amounts of the other light elements are produced: 
D and '^He at the level of about 10"^ by number, 
and ^Li at the level of 10"^" by number. 
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Figure 1 . The light element abundances from big 
bang nucleosynthesis as a function of 7710 = lO^'^rj. 

The resulting abundances of the light elements 
are shown in Figure 1. The curves for the ^He 
mass fraction, Y, bracket the computed range 
based on the uncertainty of the neutron mean- 
life which has been taken as |^ t„ = 887 ± 2 
s. Uncertainties in the produced ''Li abundances 
have been adopted from the results in Hata et al. 
S. Uncertainties in D and ^Hc production arc 



negligible on the scale of this figure. The boxes 
correspond to the observed abundances and will 
be discussed below. 

2. Abundances 
2.1. ^He 

There are over 70 different low metallicity ex- 
tragalactic HH (ionized hydrogen) regions with 
abundance information on the "^He mass fraction, 
Y, 0/H, and N/H |,|§. The observation of 
the heavy elements is important as the helium 
mass fraction observed in these HH regions has 
been augmented by some stellar processing, the 
degree to which is given by the oxygen and nitro- 
gen abundances. In an extensive study based on 
the data, it was found that the data is well rep- 
resented by a linear correlation for y vs. O /H and 
Y vs. N/H. It is then expected that the primor- 
dial abundance of ^He can be determined from 
the intercept of that relation. The '*He mass frac- 
tion [H based on 62 distinct HII regions is 

Fp = 0.234 ±0.002 ±0.005 (1) 

The first uncertainty is purely statistical and the 
second uncertainty is an estimate of the system- 
atic uncertainty in the primordial abundance de- 
termination. The solid box for ^He in Figure 1 
represents the range (at 2CTstat) from (|^). The 
dashed box extends this by including agys- A 
somewhat lower primordial abundance of Yp = 
0.230 ± .003 is found by restricting to the 32 most 
metal poor regions A similar result is found 
using a Bayesian analysis which does not depend 
on a linear regression of the data Q . 

There may be some evidence however, that the 
^He mass fraction is somewhat higher than in Eq. 
(|l|). First, there is a potential problem concern- 
ing underlying stellar absorption [p^ with one 
data point (I Zw 18 NW) that is at low 0/H 
and has small errors due its having been observed 
more than once. When this point is removed, 
Yp becomes 0.237 ± 0.003 for 61 data points. 
There is also a new (preliminary) compilation of 
data from 45 HII regions by Izotov and Thuan 
||ll| which includes new observations as well as 
the ones published in Based on their "self- 
consistent" analysis of 5 helium lines, they find 
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Yp = 0.244 ± 0.002, a significantly higher value 
for Yp. However, some of the lines are very sen- 
sitive to the electron density in the HII region. 
When densities determined from sulphur obser- 
vations are used (as in they find a lower 
(but still high) value of Yp = 0.239 ± .002. This 
is to be compared with Yp = 0.0234 ± 0.003 for 
the data in from ||,| with I Zw 18 NW re- 
moved. When all of the available data is used, 
I find Yp = 0.238 ± 0.002 though this should be 
taken as a tentative value. I will comment below 
on the effect of higher Yp on the BBN consistency 
analysis. 

2.2. ^Li 

The ^Li abundance is also reasonably well 
known. In old, hot, population-II stars, ^Li is 
found to have a very nearly uniform abundance. 
For stars with a surface temperature T > 5500 K 
and a metallicity less than about l/20th solar (so 
that effects such as stellar convection may not be 
important), the abundances show little or no dis- 
persion beyond that which is consistent with the 
errors of individual measurements. There is ^Li 
data from nearly 100 halo stars, from a variety of 
sources. I will use the value given in ||l^ as the 
best estimate for the mean ^Li abundance and its 
statistical uncertainty in halo stars 



Li/H : 



(1.6±0.ll°;™)xl0- 
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(2) 



The first error is statistical, and the second is 
a systematic uncertainty that covers the range 
of abundances derived by various methods and 
is probably an overestimate given the analysis in 
12 1 . The solid box for ^Li in Figure 1 represents 
the 2a range including this uncertainty. The third 
set of errors in Eq. (||) accounts for the possibil- 
ity that as much as half of the primordial ^Li has 
been destroyed in stars, and that as much as 30% 
of the observed ^Li may have been produced in 
cosmic ray collisions rather than in the big bang. 
The dashed box in Figure 1 , accounts for this ad- 
ditional uncertainty. Observations of ^Li, Be, and 
B help constrain the degree to which these effects 
play a role |Q. For ^Li, the uncertainties are 
clearly dominated by systematic effects. 



2.3. D and ^He 

We have three basic types of abundance in- 
formation on these two light elements: 1) ISM 
data, 2) solar system information, and perhaps 
3) a primordial abundance of D/H from quasar 
absorption systems. The best measurement for 
ISM D/H is (l| 

(D/H)isM = 1.60 ± Omt° fo X 10"^ (3) 



The lower bound using this deuterium value 
establishes an upper bound on r/ which is 
robust and is shown by the lower right of 
the solid box in Figure 1. The solar 
abundance of D/H is inferred from measure- 
ments of ^He in meteorites and the solar 
wind. These measurements indicate that |l5| 
|6| ( ^He/H)^ = (1.5 ±0.2 ±0.3) X 10-5 and 

(D/H)0 ss (2.6 ± 0.6 ± 1.4) x 10"^. 

There have been several recent reported mea- 
surements of D/H is high redshift quasar absorp- 
tion systems. Such measurements are in princi- 
ple capable of determining the primordial value 
for D/H and hence t], because of the strong and 
monotonic dependence of D/H on r/. Though 
there are several measurements of D/H in such 
systems, a unique value has yet to be estab- 
lished. The first of these measurements ||l^ as 
well as a more recent (relatively low z) observa- 
tion indicate a rather high D/H ratio, D/H ss 
2 xlO^"'. Other high D/H ratios were reported 
in ||l^. However, there are reported low values 
of D/H in other such systems Q with values 
D/H ~ 2.5 X 10~5, significantly lower than the 
ones quoted above. (The revised value for these 
is about 3.5 xlO^^.) The range of quasar ab- 
sorber D /H is shown by the dashed box in Figure 
1. It is probably premature to use either of these 
values as the primordial D/H abundance in an 
analysis of big bang nucleosynthesis, but it is cer- 
tainly encouraging that future observations may 
soon yield a firm value for D/H. It is also very 
unlikely that both are primordial and indicate an 
inhomogeneity ||2l|] due to microwave background 
constraints. It is however important to note that 
there does seem to be a trend that over the his- 
tory of the Galaxy, the D/H ratio is decreasing, 
something we expect from galactic chemical evo- 
lution. The amount of deuterium astration that 
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has occurred is model dependent. 

Finally, there are in addition to the solar ■^He 
values, measurements of ■^He/H in the ISM with 
H ( ^He/H)jg = 2.ltl X 10-5 and in galac- 
tic HII regions which show a wide disper- 
sion which may be indicative of pollution or a 
bias H ( 3He/H)jjjj ~ 1 - 5 x 10"^ There are 
also observations of "^He in planetary nebulae ||25|| 
which show a very high '^He abundance of ■^Hc/H 
~ IQ-'^, corresponding to the production of '^He 
in low mass stars as is expected from standard 
stellar models [ p6| . 

3. Likelihood Analysis 

Monte Carlo techniques have proven to be a 
useful form of analysis for big bang nucleosyn- 
thesis p7| , p8P ] and an analysis of this sort was 
performed ]29[ concentrating on "^He and ^Li. It 
should be noted that in principle, two elements 
should be sufhcient for not only constraining the 
one parameter (rj) theory of BBN, but also for 
testing for consistency. Here, I will not describe 
the procedure for calculating the various likeli- 
hood functions but simply show the results. For 
more information, the reader is referred to [E9l. 
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Figure 2. Likelihood distribution for each of '^He 
and ^Li, shown as a function of rj. 



What are plotted in Figure || are the individ- 
ual likelihood functions for '^He and ''Li as a func- 
tion of 7710. The assumed observational element 
abundances were taken from Eqs. (§) and (|). 



The one-peak structure of the "^He curve corre- 
sponds to its monotonic increase with 77, while the 
two-peaks for ''Li arise from its passing through 
a minimum ( i.e., a given observed value of ''Li, 
there are two likely values of 77) as can be surmised 
from Figure 1. As one can see there is very good 
agreement between ^He and ^Li in the vicinity of 
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Figure 3. Combined likelihood for simultaneously 
fitting ''He and ^Li, as a function of 77. 

The combined likelihood, for fitting both ele- 
ments simultaneously, is given by the product of 
the two functions in Figure ^ and is shown in Fig- 
ure ||. The likelihood prediction (95% CL) for r] 
is 7710 = 1.8^^4* corresponding to 



.0061;°°? 



(4) 



where VIb is the fraction of critical density in 
baryons, and h is the Hubble parameter scaled 
to 100 km/s/Mpc. 

It is interesting to note that the central (and 
strongly) peaked value of 7710 determined from the 
combined ^He and^Li likelihoods corresponds to 
avalueof D/H ~ 1.8 xlO""*, very close §3] to the 
high value of D/H in quasar absorbers |17, p"8|JlSt| . 
This is one of the main motivations for studying 
galactic chemical evolution models with high ini- 
tial D/H and strong D destruction histories. 

The effects of the higher values of Yp discussed 
in section 2.1 can be seen in Figure ^. Shown as 
solid curves are the same two likelihood functions 
from Figure ^ for "^He and ''Li. Also shown in solid 
are the likelihood functions for high D/H (at low 
7710) and low D/H (at high 7710). As can be seen, 



5 



there is striking agreement between *He, ^Li, and 
high D/H. The dotted curve corresponds to Yp ~ 
0.238±0.002 and as one can see there is still broad 
agreement with ^Li and high D/H. The maximum 
likelihood value for ?/ is 7710 = 1.9 with a 95% 
CL range of 1.6 < 7710 < 4.6. The dashed curve 
shows the result for Yp = 0.244 ± 0.002. In this 
case, we learn very little from ''He as its likelihood 
function is now spread over nearly all values of 77 
in the region of interest. There is overlap between 
this distribution and ^Li (though it now favors 
the high rj peak) as well as overlap with both high 
and low D/H. The maximum likelihood value for 
77 in this case is 7710 — 3.7 with a 95% CL range 
of 1.8 < 7710 < 5.0 from both "'He and '^Li. For 
completeness, I note that the 95% CL range for 
high D/H is 1.3 < 7710 < 3.8 and for low D/H it 
is 4.8 < 7yio < 6.2. 
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Figure 4. A comparison of likelihood distribu- 
tions for "He with Yp = 0.234, 0.238, and 0.244 
and "^Li, as well as high and low D/H. 



4. Implications 

4.1. Chemical Evolution 

Because we can not directly measure the pri- 
mordial abundances of any of the light element 
isotopes, we are required to make some assump- 
tions concerning the evolution of these isotopes. 
As has been discussed above, ''He is produced in 
stars along with oxygen and nitrogen. ^Li can 
be destroyed in stars and produced in several 



(though still uncertain) environments. D is to- 
tally destroyed in the star formation process and 
•^He is both produced and destroyed in stars with 
fairly uncertain yields. It is therefore preferable, 
if possible to observe the light element isotopes in 
a low metallicity environment. Such is the case 
with *He and ^Li, and we can be fairly assured 
that the abundance determinations of these iso- 
topes are close to primordial. If the quasar ab- 
sorption system measurements of D/H stabilize, 
then this too may be very close to a primordial 
measurement. Otherwise, to match the solar and 
present abundances of D and '^He to their primor- 
dial values requires a model of galactic chemical 
evolution. 

It is beyond the scope of this discussion to re- 
view galactic chemical evolution, but some com- 
ments are in order. Deuterium is always a mono- 
tonically decreasing function of time in chemical 
evolution models. The degree to which D is de- 
stroyed, is however a model dependent question. 
The evolution of '^He is considerably more com- 
plicated. Stellar models predict that substantial 
amounts of '^He are produced in stars between 1 
and 3 Mq and is consistent with the observation 
of high •^He/H in planetary nebulae How- 
ever, the inclusion of the standard '^He yield in 
chemical evolution models leads to an overpro- 
duction of ^He/H particularly at the solar epoch 
||2j,^. This problem is compounded in mod- 
els with an intense period of D destruction. In 
Scully et al. |Q , a dynamically generated super- 
novae wind model was coupled to models of galac- 
tic chemical evolution with the aim of reducing 
a primordial D/H abundance of 2 xlO^'' to the 
present ISM value without overproducing heavy 
elements and remaining within the other observa- 
tional constraints typically imposed on such mod- 
els. The time evolution of D/H and '^He/H from 
this model is shown in Figure 5. The dashed curve 
assumes standard yields for '^He. 

The overproduction of '^He relative to the solar 
metcoritic value seems to be a generic feature of 
chemical evolution models when "^He production 
in low mass stars is included. This conclusion is 
independent of the chemical evolution model and 
is directly related to the assumed stellar yields of 
•^He. It has recently been suggested that at least 
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some low mass stars may indeed be net destroy- 
ers of ^He if one includes the effects of rotational 
mixing in low mass stars on the red giant branch 
@,||. It was shown H that these reduced ^He 
yields in low mass stars can account for the rel- 
atively low solar and present day ^He/H abun- 
dances observed as shown by the solid '^He curve 
in Figure 5. 




Figure 5. The evolution of D/H and ■^He/H as a 
function of time. 



4.2. Dark Matter 

The implications of the resulting predictions 
from big bang nucleosynthesis on dark matter are 
clear. First, if = 1 (as predicted by inflation), 
and Qb ^0.1 which is certainly a robust conclu- 
sion based on D/H, then non-baryonic dark mat- 
ter is a necessity. Second, on the scale of small 
groups of galaxies, > 0.05, and is expected to 
sample the dark matter in galactic halos. This 
value is probably larger than the best estimate 
for from equation (^. VlBh? = 0.006 corre- 
sponds to riB — 0.025 for h = 1/2. In this event, 
some non-baryonic dark matter in galactic halos 
is required. This conclusion is unchanged by the 
inclusion of the high D/H measurements in QSO 
absorbers. In contrast, the low D/H measure- 
ments would imply that ils/i^ = 0.023 allowing 
for the possibility that ^Ib — 0.1. In this case, no 
non-baryonic dark matter is required in galactic 
halos. 

4.3. Particle Properties 

Limits on particle physics beyond the stan- 
dard model from BBN are mostly sensitive to the 
bounds imposed on the *He abundance. As is well 



known, the '*He abundance is predominantly de- 
termined by the neutron-to-proton ratio just prior 
to nucleosynthesis and is easily estimated assum- 
ing that all neutrons are incorporated into ^He. 
As discussed earlier, the neutron-to-proton ratio 
is fixed by its equilibrium value at the freeze-out 
of the weak interaction rates at a temperature 
jy 1 MeV modulo the occasional free neutron 
decay. Freeze-out is determined by the competi- 
tion between the weak interaction rates and the 
expansion rate of the Universe 



^T/ ^ r„k(T/) = HiTf) ^ ^GNNTf^ (5) 

where N counts the total (equivalent) number 
of relativistic particle species. The limit on 
comes about via the change in the expansion rate 
given by the Hubble parameter. 



■p = [Asm 



3 " 90 ■ 8"' ''^^ 

when compared to the weak interaction rates. 
Here A"sm refers to the standard model value for 
N. At T - 1 MeV, A^sm = 43/4. The presence 
of additional neutrino flavors (or any other rela- 
tivistic species) at the time of nucleosynthesis in- 
creases the overall energy density of the Universe 
and hence the expansion rate leading to a larger 
value of Tf, (n/p), and ultimately Yp. Because 
of the form of Eq. ^ it is clear that similar to 
placing limits | |3^ on A'^, any changes in the weak 
or gravitational coupling constants can be also 
constrained (for a recent discussion see Q). 

Just as ''He and ^Li were sufficient to deter- 
mine a value for rj, a limit on N^, can be obtained 
as well Since the light element abun- 

dances can be computed as functions of both ij 
and Ni,, the likelihood method can be extended 
to include both rj and Ni, as arguments. A three- 
dimensional view of the combined likelihood func- 
tions is shown in Figure |^. In this case the 
high and low rj maxima of Figure H, show up as 
peaks in the L — rj — space. The peaks of the 
distribution as well as the allowed ranges of rj and 
N„ are more easily discerned in the contour plot 
of Figure @ which shows the 50%, 68% and 95% 
confidence level contours in the likelihood func- 
tion. The crosses show the location of the peaks of 
the likelihood functions. L47 peaks at A"^ = 3.0, 
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?7io = 1-8 (in agreement with our previous results 
(^) and at = 2.3, 7710 = 3.6. The 95% con- 
fidence level allows the following ranges in 77 and 

1.6 <N^< 4.0 1.3 < 7710 < 5.0 (7) 




account the reduced abundance of these parti- 
cles at the time of nuclcosynthesis|^. For ex- 
ample, when the limit is applied to scalar fields 
we find that a single additional scalar degree of 
freedom (which counts as |) such as a majoron is 
allowed. On the other hand, in models with right- 
handed interactions, and three right-handed neu- 
trinos, the constraint is severe. The right-handed 
states must have decoupled early enough to en- 
sure that their present temperature have the ratio 
{T„j^/T^^y < 1/3. If right-handed interactions 
are mediated by additional gauge interactions, as- 
sociated with some scale Mz' , and if we assume 
that the right handed interactions scale as M^, 
with a standard model-like coupling, then it is 
possible to set a limit on the mass scale Mz' ^1.6 
TeV! In general this constraint is very sensitive to 
the BBN limit on 7V^. Of course BBN can be used 
to constrain much more than the number of neu- 
trino flavors and the reader is referred to |Q for 
a more complete review. 



Figure 6. The combined two-dimensional likeli- 
hood function, L47 for simultaneously fitting ''He 
and ''Li as functions of both rj and N^. 

Since it is the '^Hc abundance which is most 
sensitive to iVV, one might expect that the limit 
on depends on our choice of 1^. It turns out 
however, that the upper limit to N^, is relatively 
firm. For Yp = 0.238 ± 0.002, the peak value of 
rjio = 1.9 corresponds to a peak value N^, = 3.2 
and a 95% CL upper limit, < 4.1. Similarly, 
for Yp = 0.244±0.002, the peak value of 7710 = 3.7 
corresponds to a peak value N^, — 3.0 and a 95% 
CL upper limit, iV^ < 4.0. Although the value of 
Yp is moving up allowing for more neutrinos, the 
peak value of 77 is also moving up and compen- 
sates by allowing fewer neutrino fiavors. The net 
result is a robust limit. 

One should recall that the limit derived above 
is not meant for neutrinos in the strictest sense. 
That is, the limit is only useful when applied 
to additional particle degrees of freedom which 
necessarily do not couple to the Z°. For very 
weakly interacting particles, one must take into 




Figure 7. Likelihood contours in the combined 
likelihood function for ^He and ^Li. The con- 
tours represent 50% (innermost), 68% and 95% 
(outermost) confidence levels. The crosses mark 
the points of maximum likelihood. 
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